Non-analytical Angular Dependence of the Upper Critical Magnetic Field in a 

Quasi-One-Dimensional Superconductor 
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We have derived the so-called gap equation, which determines the upper critical magnetic field, 
perpendicular to conducting chains of a quasi-one-dimensional superconductor. By analyzing this 
equation at low temperatures, we have found that the calculated angular dependence of the upper 
critical magnetic field is qualitatively different than that in the so-called effective mass model. In 
particular, our theory predicts a non-analytical angular dependence of the upper critical magnetic 
field, H c2 (0) - H c2 (a) ~ a 3/2 , when magnetic field is close to some special crystallographic axis and 
makes an angle a with it. We discuss possible experiments on the superconductor (DMET)2l3 to 
discover this non-analytical dependence. 

PACS numbers: 74.70.Kn, 74.25.Op 
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Upper critical magnetic field, which corresponds to de- 
struction of superconductivity in type II superconduc- 
tors, is known to be one of the most fundamental prop- 
erties of a superconducting state. The first calculations of 
the upper critical magnetic field were done in the frame- 
work of the phenomenological Ginzburg-Landau (GL) 
theory (see, for example, [1,2]) before the creation of the 
Bardeen-Cooper-Schrieffer (BCS) microscopic theory of 
superconductivity. Later, it was shown [3] that the GL 
theory is a limiting case of the BCS theory at T c — T <C T c 
and the upper critical magnetic fields were calculated [4] 
at T c — T -C T c for superconductors with anisotropic 
electron spectra, where T c is superconducting transition 
temperature in the absence of a magnetic field. Using the 
microscopic Gor'kov equations, the upper critical field 
was calculated for a 3D isotropic superconductor at zero 
temperature [5] and at arbitrary temperatures [6] . As for 
superconductors with anisotropic electron spectra, the 
common belief is that we can apply the results [4], ob- 
tained at T c — T <C T c , at any temperature, including 
T <C T c . The results [4] are usually called the effective 
mass (EM) model. 

The main goal of our Letter is to show that the shape 
and topology of the Fermi surface (FS) play a crucial 
role in determination of angular dependence of the upper 
critical magnetic field at low temperatures. We consider 
a quasi-one-dimensional (Q1D) superconductor, which is 
characterized by two open slightly corrugated sheets of 
the FS. By using the Gor'kov equations [3], we derive 
the so-called gap equation, determining the upper crit- 
ical magnetic field, perpendicular to conducting chains 
in a Q1D superconductor. As a result, we obtain a 
rather complicated integral equation, which we numer- 
ically solve at T -C T c . Our numerical analysis of this 
integral equation shows that the EM model cannot be 
applied to Q1D case at T <C T c even at qualitative level. 
Our main finding is that we predict non-analytical an- 
gular dependence of the upper critical magnetic field, 
H C 2(0) — H C 2(a) ~ a 3 / 2 , in the case, where magnetic field 
is close to some special crystallographic axis and makes 



an angle a with it. This fact is in a sharp disagreement 
with a common belief, based on the results of the EM 
model, that H C 2(0) — H c2 (a) has to be proportional to a 2 . 
Our second finding is that superconducting nuclei (i.e., 
solutions of the gap integral equation) are not of an expo- 
nential shape. We show that they decay very slowly and 
change their signs with distance. It is important that the 
above described phenomena are novel and due to quasi- 
classical effects of an electron motion in a magnetic field 
along open sheets of the Q1D FS in a single Brillouin 
zone. They are different from quantum effects of an elec- 
tron motion in the extended Brillouin zone, considered 
in Refs.[7,8]. Moreover, for discovery of non-analytical 
angular dependence, we need different experimental con- 
ditions than for investigation of the so-called Reentrant 
superconductivity [7-11]. We propose to investigate ef- 
fects, suggested in the Letter, in the Q1D superconduc- 
tor (DMET)2l3, where the upper critical magnetic fields 
have been recently measured along all three principal di- 
rections [12]. It has also been pointed out [12] that super- 
conductivity in the above mentioned compound is very 
far from the Reentrant superconducting regime [7], in 
contrast to superconductivity in (TMTSF^X materials 
[7-11]. 

Let us consider a superconductor with the following 
Q1D electron spectrum, 

Se ± {p) = ±v F (p x TPf) - 2t y cos(p y a y ) ~ 2t z cos(p z a z ), 

(1) 

in a magnetic field, 



H = (0, H cos a, H sin a) , A = (0, Hx sin a, — Hx cos a) , 

(2) 

perpendicular to its conducting chains. [Here, +(-) 

j ty ^* t z 

.y oiiu a z crystallo- 
graphic axes; vf andp^ are the Fermi velocity and Fermi 
momentum, respectively; h = 1.] 

To determine electron wave functions in the mixed rep- 



stands for right (left) sheet of the Q1D FS (1;, w y 
are electron hoping integrals along a y and a z 
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resentation, ^f(x,p y ,p z ), where 

^t(x,y,z) = exp[ip F x}exp(ip y y)exp(ip z z)^f(x,p y ,p z ), 

(3) 

we use the so-called Peierls substitution method, p x =p 
p F -> -id/dx, p y ^ p y - eAy/c, p z -> p z - eA z /c. As 
a result, we obtain the following electron Hamiltonian in 
the presence of a magnetic field: 



// | "'i ~j~ 2ty cos | Pit' 1 It 



vf' 

oj z 



- 2t z cos (p z a z H x 



(4) 



where u; y = evpHdySma/c and w z = evpHa z cos a/c. 

In this Letter, we ignore quantum effects of an electron 
motion in a magnetic field in the extended Brillouin zone 
[7-11] and use the so-called eikonal approximation [3]. 
Note that we consider the case of small angles, a <C 1, 
where oo z >• cj y , which is important for non-analytical de- 
pendence of the upper critical field. As shown in Ref.[7], 
the quantum effects are small only at high enough tem- 
perature, where 



T> T* 



u z 
8tt 2 



(5) 



(see Eq.(6) of Rcf. [7]). Under condition (5), we can lin- 
earize the Hamiltonian (4) with respect to a magnetic 
field, 

H = TivF-r — 2t y cos(p y a y ) — 2t y s'm(p y a y ) 
ax v f 

LO z X 

- 2t z cos(p z a z ) + 2t z sin(p z a z ).(6) 

vf 

It is important that the corresponding Schrodinger equa- 
tion for wave functions in the mixed representation, 



H*S> e {x,Py,Pz) = 5et> e (x,p y ,p z ), 
can be exactly solved: 



(7) 



ty € (x,p y ,p z ) — exp(±iSex/v F ) cxp[±i(j> y {p y , x)] 

exp[±i(j) z (p z ,x)}, (8) 



where 



<t>y(Py, x) = — cos(p y a y )x + sin(p y a y )x 2 , 
vf v f 



j f 

4>z(Pz,x) = — cos(p z a z )x- -4^- sin(p z a z )x 2 . (9) 

VF V F 

Since the electron spectrum and wave functions are 
known, the corresponding finite temperatures Green 
functions can be determined by means of the standard 
procedure [13]: 



V F 



x exp[ip y (y - yi)] exp[ip z (z - zi)] exp 

x exp[±i2i y cos(p y a y )(x — xi)/vp] 
x exp[±i2i z cos(p z a z )(x — x\)/vf] 
x cxp[±it y Lu y sm(p y a y )(x 2 — x 2 )/v F ] 
x exp[=Fii z w z sin(p z a z )(x 2 - x\)/v 



Tu n (x - Xi) 
V F 



F\- 



(10) 



The so-called gap equation, determining superconduct- 
ing transition temperature in the presence of the mag- 
netic field (2), can be derived by using the Gor'kov equa- 
tions for non-uniform superconductivity [14]. As a result, 
we obtain: 



A(») = | 



2irTdx\ 



Jo 



\x— x\ \ >d 

(x^-xl) 



vf sinh(27rT|a: — x\\/vf) 



2t z LO z 



(x^xl) 



A( Xl )(ll) 



where g is an effective electron coupling constant, d is 
a cutoff distance [15]. Here, we rewrite Eq.(ll) in more 
convenient way: 



Jo 



2irTdz 



\z\>d V F sinh(27rT|z|/w F ) 
2t z u 



Jo 



z(z + 2x) 



A(x + Q2) 



[Note that the Pauli paramagnetic spin-splitting effects 
are ignored in all equations above, which means that 
the upper critical magnetic field is supposed to be 
much smaller than the so-called Clogston-Chandrasekhar 
paramagnetic limit [17]. Such situation, for example, 
is experimentally realized in the Q1D superconductor 
(DMET) 2 I 3 [12].] 

Let us determine the GL slope of the upper critical 
magnetic field in the vicinity of superconducting tran- 
sition temperature. To achieve this goal, we need to 
take into account that in the GL region, (T c -T)«T C , 
vf/2itT c <C vp/ \/t y LiJy, vf / ^Jt z ui z . In this case, we can 
expand the integral equation (12) in terms of small pa- 
rameter, \z\ ~ vf/2ttT c . As a result of such expansion 
procedure, we obtain the following differential equation: 



d 2 A(x) 2 8(t 2 A| + i|w|) 



dx 2 



A(x) 



f 

Jo 



wT r z 2 dz 



vf sinh(27rT c /?;ir) 



+ 



1 



2nTdz 



i d vf smh(2nT z / v f) 
If we take into account that 

2irT c dz 



A(x) = 0. 



(13) 
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Vf sinh(27r T c z/vf) ' 
then we can rewrite Eq.(13) in the following way: 



(14) 



2 d 2 A(x) t (2ttH\ 



Py,Pz 



-tA{x) = 0, (15) 
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where 



f 2 

S X 



[Here 



7C(3)4 
16(ttT c ) 2 



£ 2 



7C(3)t 



2 a 2 



8(7tT c 



7C(3)go| 
8(7tT c ) 2 

T c -T 



(16) 



= irhc/e is the flux quantum, £ x , £ y , and £ z are 
the coherence lengths along a^, a y , and a 2 axes, corre- 



spondingly.] 
tionship: 



Note that above we use the following rela- 



z dz 
sinh(z) 



7C(3) 




(17) 



where C( n ) is the Rcimann zeta function [18]. To find 
the GL slope of the upper critical magnetic field, per- 
pendicular the conducting chains, we need to determine 
the lowest energy level of the Schrodinger-like GL equa- 
tion (15). As a result, we obtain 



FIG. 1: A typical solution of Eqs.(20),(21) is shown for a = 
0. Note that it is an oscillatory slow decaying function of 
coordinate x, in contrast to exponential solution of Eq.(15) in 
the EM model [4]. 



1 



sin 2 a 



HMO)' 



(18) 



Hc 2 (a)-Hc 2 (0) 



-0.05 



where 



H c2 (0,T) 



T c -T 
T c 

T c -T 



Tc — T 



4V2^ 2 cT 2 
7((3)ev F t z a z \ T c 
4V2tt 2 cT 2 (Tp- 



-0.10 



-0.15 



7£(3)evFt y a y \ T r 




a (degrees) 



[Note that Eq.(18) is usually called EM model and ap- 
plied to fit the experimental upper critical magnetic fields 
at any temperature, including T -C T c . On the other 
hand, we pay attention that Eqs.(13)-(19) are derived 
under the GL condition (T — T c ) <C T c , which is equiv- 
alent to the following two conditions: T 2 t z u z and 
T 2 3> tyUjy. It is important that the latter inequal- 
ities can be rewritten as: H -c H C 2(0,T — 0) and 
Hsma < H c2 (ir/2,T = 0).] 

Below, we consider the gap equation (12) at low tem- 
perature, T c 3> T T* , where we can formally employ 
T = in Eq.(12): 



FIG. 2: Angular dependence of the upper critical mag- 
netic field at T = 0, H C 2{a) — H C 2(0), calculated from 
Eqs.(20),(21), is shown by dots, which are well fitted by func- 
tion -a 3/2 (solid line). The EM model result (18), where 



H c2 {a) -H c2 (0) 
parison. 



-a , is shown by dashed line for a com- 



A(x) 



fJ 



dz 



Jo[/3 sin(a)z(2x + z)] 



2 J y>2LQ t z d/v F Z 

x Jo[cos(cv)z(2a; + z)]A(x + z) + Jq[/3 sm(a)z(2x — z)] 
xj [cos{a)z{2x- z)]A(x- z)\, (20) 

where j3 — t y a y /t z a z , ujq — evpHa z /c. It is important 
that the effective electron coupling constant, <?, and cutoff 
distance, d, can be eliminated from Eq.(20) by using the 
following relationship: 



1 



dz 



nT c d/v F v F sinh(z) : 



(21) 



which is a result of Eq.(14). 

Note that experimental value of the parameter j3 in 
(DMET)2l3 superconductor is estimated as /? ~ 10 [12]. 
Below, we analyze Eqs.(20),(21) numerically by solving 
the gap integral equation (20) under the condition (21) 
for f3 = 10. Let us first consider the case a = 0, where 
magnetic field is applied along & y axis. A typical solution 
of Eq.(20), which is called superconducting nucleus, in 
this case is shown in Fig.l. As seen from Fig.l, in our case 
superconducting nucleus changes its sign and slowly de- 
cays in space, in contrast to the results of the EM model 
[4]. Note that ata^O solutions of Eqs.(20),(21) become 
more complicated, but they retain the above mentioned 
unusual properties. In Fig. 2, we show the calculated an- 
gular dependence H c2 (a) — H c2 (0) and its fit by function 
—a 3 / 2 . Note that the agreement between the calculated 
angular dependence and function —a 3 / 2 is very good. For 
comparison, we also show dependence (18), expected in 



the EM model, where H c2 {a) - H c2 (0) 



[19]. 
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Hc 2 (a)-Hc 2 (0) 
H EM c 2 (a)-H EM c 2 (0) 




a (degrees) 



FIG. 3: Calculated by means of Eqs.(20),(21) and normalized 
angular dependence of the upper critical magnetic field (see 
the main text). 



Hc 2 (a)-H^(a) 




FIG. 4: Calculated difference of angular dependence of the 
upper critical magnetic field, given by Eqs.(20),(21), and that, 
given by the EM model (18), H c2 (a) - H? 2 M (a). 



In Fig. 3, we plot the calculated angular dependence 
of the upper critical magnetic field, normalized on the 
corresponding result (18) of the EM model, [H C 2(a) — 
H c2 {Q)]/[Hg M (a) - Hg M (Q)]. As it follows from Fig.3, 
the maximum deviations from the EM model occur at 
low angles and in the vicinity of some angle a ~ 5°. 
At low angles, the calculated in the Letter upper critical 
magnetic field exhibits different angular dependence than 
that in the EM model (18), as discussed above. To clarify 
nature of the minimum in Fig.3 at a ~ 5°, we plot the 
difference, H C 2(a) — H^ 2 M (ct), in Fig. 4. As seen from 
Fig. 4, the maximum difference corresponds to a ~ 5.6° - 
angle, which we relate to the following theoretical value: 

a* = arctan(l//3) = arctan(l/10) ~ 5.71°. (22) 

Note that, under condition (22), both Bessel functions 
in Eq.(20) have the same arguments and some kind of 
resonance appears. We suggest to measure experimen- 
tally the position of the peak in the angular depen- 
dence H C 2 (a) — Hf 2 M (a) to carefully determine the ratio 
P = t y a y /t z a z from Eq. (22). 

To summarize, we have shown that the EM model [4] is 
not adequate to describe the upper critical magnetic field 
in superconductors with anisotropic electron spectra at 
low temperatures. For the case of a Q1D superconductor, 



we have found non-analytical angular dependence of the 
upper critical magnetic field, H C 2(a) — H C 2(0) ~ — a 3 / 2 , 
where a magnetic field is perpendicular to conducting 
axis, sl x , and makes angle a with axis & y . In addi- 
tion, some angular resonance is predicted for "magic" 
direction of a magnetic field (22). We suggest to test 
the above mentioned predictions of the Letter on the 
Q1D superconductor (DMET)2l3, where the upper crit- 
ical magnetic fields along the main crystallographic axes 
have been recently measured [12]. In our opinion, un- 
conventional shapes of superconducting nuclei as well as 
the non-analytical angular behavior of the upper critical 
field, found in the Letter, may reflect the existence of un- 
usual vortex lattice in Q1D superconductors. Therefore, 
we also suggest experimental studies of the vortex lattice 
at magnetic fields, corresponding to small values of angle 
a in Eq.(2). 

Let us prove that the (DMET)2l3 superconductor sat- 
isfies the condition of a validity of our theory, 

T c > T > T* , (23) 

at experimentally used lowest temperature, T ~ 0.05 K , 
where T* is given by Eq.(5) and T c = 0.5 K [12]. If 
we take from Rcf.[12] the typical experimental values, 
H V C 2 = °- 2 T > v f - 0.4 x I0 7 cm/s, a z = 15.8 A, t z ~lK, 
we obtain T* ~ 0.006 K. Therefore, we conclude that the 
suggested in the Letter theory is applicable to the super- 
conductor (DMET)2l3 at the lowest experimental tem- 
perature [12]. Note that, for neglecting the quantum cor- 
rections [7,8] and, thus, the Reentrant Superconductivity 
effects [7-11], it is also important that 4i z /w z ~ 27 ^> 1, 
as has been already mentioned in Rcf.[12]. 

We point out that in a geometry, considered in the Let- 
ter, experiments were performed in Ref . [20] in the super- 
conductor (TMTSF) 2 C10 4 in low magnetic fields, H < 
H C 2 , to demonstrate another phenomenon - the so-called 
lock-in effect. To avoid lock-in effect [20], the experi- 
ments, suggested by us, have to be performed at magnetic 
fields, which satisfy the condition H z ~ H y c2 sin a ^> 
[20]. Although is not known in the superconductor 
(DMET)2l3, it is clear that in this typical type-II super- 
conductor iJJj -C H* 2 — 0.02 T, which shows that the 
above mentioned condition is presumably satisfied. 
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